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It is said that the microgravity environment positively affects the quality of 
protein crystal growth. The formation of a protein depletion zone and an 
impurity depletion zone due to the suppression of convection flow were thought 
to be the major reasons. In microgravity, the incorporation of molecules into a 
crystal largely depends on diffusive transport, so the incorporated molecules will 
be allocated in an orderly manner and the impurity uptake will be suppressed, 
resulting in highly ordered crystals. Previously, these effects were numerically 
studied in a steady state using a simplified model and it was determined that the 
combination of the diffusion coefficient of the protein molecule (D) and the 
kinetic constant for the protein molecule (P) could be used as an index of the 
extent of these depletion zones. In this report, numerical analysis of these 
depletion zones around a growing crystal in a non-steady {i.e. transient) state is 
introduced, suggesting that this model may be used for the quantitative analysis 
of these depletion zones in the microgravity environment. 



Keywords: microgravity; transient and homogeneous field; transient and diffusive field; 
protein crystal; numerical model; protein depletion zone; impurity depletion zone. 



1. Introduction 

Protein crystal growth experiments are a promising area in 
the usage of microgravity to contribute to structural biology 
(McPherson, 1999; Littke & John, 1986; Kundrot et al, 2001; 
Vergara et a/., 2003). When protein molecules are taken into a 
crystal, a spherical area of low protein concentration is formed 
around the growing crystal. In the terrestrial environment, a 
density-driven flow occurs to supply protein molecules to the 
low concentration area, disturbing this area. However, in 
microgravity, this density-driven flow does not occur, so the 
protein molecules are supplied to the crystal only by thermal 
diffusion caused by Brownian motion. Therefore, the low 
protein concentration area around the growing crystal is 
maintained, resulting in the formation of a protein concen- 
tration depletion zone (protein depletion zone, PDZ). The 
PDZ helps grow the crystal at a low supersaturation (Chernov, 
1998; Otalora et al, 2001), eventually suppressing the disorder 
of protein molecules in the crystal. Following similar steps as 
the protein molecules, a low-impurity concentration area 
around the growing crystal (impurity depletion zone, IDZ) is 
formed (Chernov, 1998; Thomas et al, 2000), also suppressing 
crystal disorder. 



The PDZ and IDZ in microgravity in the steady and 
diffusive states have been analyzed and discussed (Tanaka et 
al, 2004). A steady state occurs when there is a constant 
concentration of the protein in a solution far from a growing 
crystal. The diffusive state occurs when the molecules diffuse 
in a convection-free environment. We have compared the 
extent of PDZ and IDZ formation in steady and diffusive 
states with the formations in steady and homogeneous states 
and discussed how the PDZ and IDZ affected protein crystal 
growth in microgravity (Tanaka et al, 2004; Inaka et al, 2012). 

However, in reality, the process of crystal growth in a 
conventional protein crystal growth experiment occurs in a 
non-steady transient state, decreasing protein concentration in 
the solution as the protein molecules are incorporated into the 
crystal. Therefore, it would be valuable to have the ability to 
employ a mathematical model corresponding to the transient 
state to know under what conditions the crystals really grow 
and how much impurity is taken into the crystals, quantita- 
tively. 

Here we introduce a rather simple numerical calculation 
model for expressing the transient state. We can compare the 
transient and homogeneous state (terrestrial gravity) with the 
transient and diffusive state (microgravity), and propose that 
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this model can be applied to the examination of the process of 
protein crystallization both in microgravity and terrestrial 
environments, quantitatively. 

2. Calculation model 

2.1 . Model of the field around a growing crystal 

There are many kinds of crystallization methods such as the 
vapour-diffusion method, dialysis, and the counter-diffusion 
method, etc. In our calculation model, for simplicity, we 
assume the batch method is used to crystallize protein because 
in this method the precipitant concentration and protein 
solubility do not change from the beginning to the end of 
crystallization. As for the nucleation, we assume that all nuclei 
start growing simultaneously and the final size of the crystal is 
the same for all crystals. In this model, we assume that crys- 
tallization occurs in a virtual sphere of radius L as shown in 
Fig. 1 and that the shape of the crystal is a sphere with a radius 
of R to simplify the calculation. 

When the final size of the crystal is R(oo), 



Xi{t) = AnRitf fii Ci(t) dt, 



(5) 



L 3 = {(« - Ce)/[C(0) - Ce]}i?(oo) 3 , 



(1) 



where C(0) and Ce are the protein concentration at the 
beginning and at the end (solubility) of crystallization, 
respectively, and n is the weight density of the crystal. 
The velocity of crystal growth is (Chernov, 1998) 



dR(t)/dt = pco[C(t)' - Ce'], 



(2) 



where C(t)' is the number of protein molecules in a unit of 
volume (1/cm 3 ) at time t after crystallization has occurred, Ce' 
is the number of protein molecules in a unit of volume at the 
protein solubility concentration and ft is the kinetic constant 
of a protein molecule. a> is the volume for one molecule and 
can be defined as a> = M/(nN A ), where M is the protein 
molecular weight and N A is Avogadro's number. 

For experimental purposes, we usually express concentra- 
tion as weight per volume at time t, C(t). Therefore, C(f)' is 
expressed as C(f)' = C{i)NJM; (2) is 
replaced with (Tanaka et al., 2004) 

dR{t)/dt = p[C(t) - Ce]/n. (3) 

Therefore, the weight of protein 
attaching to a crystal surface in an iota 
of time is 



where jii is the kinetic constant of the impurity molecule and 
Ci(t) is the impurity concentration on the surface of the crystal 
at time t. 

2.2. Transient and homogeneous model and transient and 
diffusive model 

In the transient and homogeneous model (THM) which 
represents crystal growth in the terrestrial environment, the 
concentration of the protein in the solution around the 
growing crystal is uniform and the sum of the total amount of 
the protein in the solution and in the crystal is constant. 
Therefore the concentration of the protein solution during 
crystal growth can be expressed as 



C(t) = [C(0)L 3 - R(tfn]/[L 3 - R(t) 3 ]. 



(6) 



Substituting (6) in (3), 



dR(t)/dt = (fi/n){[C(Q) - Ce] L 3 — (n — Ce)R(t) 3 } 

/[L 3 -R(ff], (7) 

Therefore, the crystal radius can be obtained from the re- 
peated calculation of (7) using the difference equation 
(Tanaka et al., 2004), 

R(t + At) = R(t) + (P/n)({[C(0) - Ce]L 3 - (n - Ce)R(t) 3 } 

/[L 3 -R(tf])At. (8) 

Considering the total amount of the impurity in the solution 
and in the crystal is constant, the concentration of the impurity 
in the solution can be obtained by the following equation, 

Ci(t + At)4it[L 3 - R(t + Atf]/3 = 

Ci(t)4jr[L 3 - R(tf]/3 - Xi{t). (9) 



., ^Virtual full sphere 



X(t) = dV(t) n = AnR{t) dR(t) n 
= AnRitf p[C(i) - Ce]dt, 



(4) 



where V(t) and R(t) are the volume and 
the radius of crystal at time t, respec- 
tively. 

When an impurity contaminates the 
solution, and if the impurity attaches to 
the surface of the crystal at a fixed ratio 
and the reverse reaction is ignored, the 
weight of the impurity attaching to the 
crystal surface in an iota of time is 




Crystal 



Virtual sphere 

(a) (b) sectioned concentrically 

by N 

Figure 1 

The conceptual configuration of the numerical model for crystal growth, (a) In actual crystallization, 
several crystals are grown in a solution, (b) To simplify this process for the model, a crystal is 
assumed to grow spherically in a virtual sphere, of which the radius, L, is related to the amount of 
the protein uptake into the crystal [see equation (1)]. To accommodate the diffusive process, the 
area between the surface of the crystal and the virtual sphere are sectioned concentrically by N. The 
diffusive processes are considered to occur between the inner section and outer section. 
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Substituting (5) into (9), the impurity concentration in the 
solution can be obtained from the repeated calculation of the 
difference equation below, 

Ci(t + At) = [Ci(t)[L 3 - R(tf] - 3R(tf fii Ci(t)At} 

/[L 3 - R(t + At) 3 ]. (10) 

On the other hand, in the transient and diffusive model 
(TDM) which represents the microgravity environment, we 
have to consider both the diffusion process in the virtual 
sphere of radius L and the crystal growth process in the centre 
of the virtual sphere. The three-dimensional diffusion equa- 
tion is dC/dt = DV 2 C, but in the case of spherical coordinates 
the partial differential equation is 

dC(r,t)/dt = D{d 2 C(r,t)/dr 2 

+ (2/r)[dC(r,t)/dr]}. (11) 

This can be applied to the area between the surface of the 
crystal and the virtual sphere of radius L. For the outer 
boundary where r is equivalent to L, no-flux of the material is 
assumed. For the boundary conditions on the surface of the 
crystal, the process of the crystal growth is 

D[dC(r, t)/dr]\ r=m = {n- C[R(t), t]} dR(t)/dt, (12) 

where C(r, t) is the protein concentration at the position of r 
from the centre of the virtual sphere at time t. The diffusion 
equation for impurities is the same as that for protein, as 
shown below with similar boundary conditions, 

dCi(r, t)/dt = Di[d 2 Ci(r, t)/dr 2 

+ (2/r)dCi(r,t)/dr], (13) 

and the process of the impurity uptake into the crystal is 

Di[dCi(r, t)/dr]\ r=R(t) = Ci[R(t), t][0i - dR(t)/dt], (14) 

where Ci(r, t) is the impurity concentration at the position of r 
from the centre of the virtual sphere at time t. 

To solve these partial differential equations numerically, it 
is common to divide a sphere along a radius into the same 
intervals to apply difference equations. However, in the case 
of a growing crystal, the sections are adsorbed into the crystal 
one after another. We may include a conditional judgment on 
whether the sections are embedded in the crystal into the 
partial differential equation. It has been found that these 
calculations may result in intolerant errors caused by the 
discontinuity of the protein concentration on the surface of 
the crystal when the section is incorporated into the crystal. 

Therefore, we divide the sphere into TV sections with a 
variable length, Ar(t), from the surface of the crystal to the 
surface of the virtual sphere whose radius is L and examine 
the diffusion in those sections. If the crystal radius is R(t) at 
time t, 

Ar(t) = [L — R(t)]/N. (15) 

The section number i (i = 1, 2, 3, . . . , N) is placed between two 
spheres whose radii are R(t) + (i — l)Ar(r) and R(t) + iAr{t). 



The process of the calculation for one time step of the 
repeated calculation of the difference equation is the 
following: 

(i) Calculate the crystal growth: based on (3), we can 
calculate the increase of the crystal radius by the following 
equation. 

Ai?(0 = ({flC(l, f) - Ce]}/n)At, (16) 

where C(i, t) is the concentration of the protein in the ;'th 
section at time t; and the radius after one time unit At is 

R(t + At) = AR(t) + R(t). (17) 

(ii) Calculate the protein concentration of the first section 
which contacts with the growing crystal: considering (12), the 
protein concentration of this section is affected by its 
adsorption onto the crystal surface and by the diffusive mass 
transfer from the next section. For the first process, the 
following relation is conserved, since the amount of protein 
which is adsorbed on the crystal is equal to the decreased 
amount of it from the section, 

[(4jr{[i?(0 + Ar(t)f - [R(t) + AR(t)f})/3]c\l, t + At) 
= (47t{[R(t) + Ar(t)f-R(tf}/3)C(1, t) 

- (4it{[R(t) + AR(t)f - R(tf}/3)n, (18) 

where C*(l, t + At) is the concentration of the protein in the 
first section after time At with the protein adsorption occur- 
ring in the first section. 

For the second process, the following difference equation is 
derived from (11) neglecting the diffusion from the inner 
section, 

C f (l, t+ At) = [C(2, t) - C(l, t)] D At/Ar(tf + C(l, t). 

(19) 

Therefore, combining these two processes, the protein 
concentration of the first section can be expressed as 

C\l, t+At)= {[R(t) + AR(t)f - R(tf}[C(l, t) - n] 
/{[R(t) + Ar(t)f-[R(t) + AR(t)f} 
+ [C(2, 0 - C(l, t)]DAt/Ar(tf + C(l, t). 

(20) 

(iii) Calculate the diffusive protein transfer between the 
surface of the crystal to the surface (i = 2 to N) of the virtual 
sphere: the difference equation for this process is derived 
from (11), 

C\i, t + At) = [C(i -1,0- C(i, t)]D[At/Ar(tf] 
x {[r 2 -2rAr{tj\/r 2 } 
+ [C(i + l,f)- C(i, t)]DAt/Ar(tf 
+ C(i, t). (21) 

(iv) Calculate the impurity concentration of the first section 
which contacts with the growing crystal: considering (14), the 
protein concentration of this section is affected by its 
adsorption onto the crystal and by the diffusive mass transfer 
from the next section. For the first process, the following 
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relation is conserved, since the amount of the impurity which 
is adsorbed on the crystal is equal to the decrease of it from 
the section, 

(4it{[R(t) + Ar(t)f - [R(t) + AR(t)f}/3)Ci\l, t + At) 

= (47t{[R(t) + Ar(t)f - R(tf}/3)Ci(l, t) - Xi(t), (22) 

where Ci(i, t) is the concentration of the impurity in the ith 
section at time t, and Ci^(l, t + At) is the concentration of the 
impurity in the first section after one time unit At. 

For the second process, the following difference equation is 
derived from (13) neglecting the diffusion from the inner 
section, 

Ci\l, t + At) = [Ci(2, t) - Ci(l, t)]DiAt/Ar(tf 

+ Ci(l, t). (23) 

Therefore, by combining these two processes the impurity 
concentration of the first section can be expressed as 

Ci\l, t + At) = ({[R(t) + AR(t)f - R(tf}Ci(l, t) 

-3R(t) 2 piCi(l,t)At) 

/{[R(t) + Ar(t)f-[R(t) + AR(t)f} 

+ [Ci(2, t) - Ci(l, t)]DiAt/Ar(tf 

+ Ci(l,t). (24) 

(v) Calculate the diffusive impurity transfer between the 
surface of the crystal to the surface of the virtual sphere: the 
difference equation for this process is derived from (13), 

Ci\i, t+At) = [Ci(i -l,t)- Ci(l, t)] 

x Di[At / Ar^tf^r 2 - 2rAr(f)]/r 2 

+ [Ci(i + 1,0- Ci(l, t)]DiAt/Ar(t) 2 

+ Ci(l,t). (25) 

(vi) Re-sectioning the sphere: after the crystal has grown 
during one time unit At, each section moves slightly to the 
outer position. Therefore, we re-divide the sphere into N 
sections with variable length Ar(t + At), 

Ar{t + At) = [L — R(t) - AR(t)]/N. (26) 

The concentration of the new section, whose number is i, is 
expressed as C(i, t + At), and can be calculated by considering 
the conservation of the amounts of the protein and the 
impurity, 

C(i, t + At) = (c\i +l,t+ At){[R(t + At) + iAr(t + At)f 

-[R(t) + iAr(t)f} 

+ C'(i, t + At){[R(t) + iAr(t)f 

- [R(t + At) + (i - l)Ar(t + A?)] 3 }) 

/ {[R(t + At) + iAr{t + At)f 

- [R(t + At) + (i - \)Ar{t + AO] 3 } , (27) 



Table 1 

Initial parameters for the calculation of the homogeneous and the 
diffusive models. 



Crystallization 


R(oc) 


Ce 


P 

(mm h~') 


D 


n 


condition 


(mm)t 


(nigral- 1 )* 


(mm 2 h-') 


(g cm" 3 ) 


Salt§ 


0.46 


4.55 


0.34 


0.360 


0.79 



t Average of the final size of the crystal radius. X Average of the final protein 
concentration. § 20 mg ml -1 purified lysozyme, 0.7 M sodium chloride, 50 mM sodium 
acetate pH 4.5. 



Ci(i, t+At)= \ Ci\i + 1, t + At){[R(t + At) + iAr(t + At)f 

-[R(t) + iAr(t)] 3 } 

+ Ci\i, t + At){[R(t) + iAr(t)f 

- [R(t + At) + (i - 1) Ar(t + At)] 3 }) 

/ {[R(t + At) + iAr(t + AO] 3 

- [R(t + At) + (i - l)Ar(f + AO] 3 }. (28) 

The calculation program was created in Microsoft C++ and 
was executed in Windows 7 or Windows XP. N was defined as 
240 and At as 2 x 1CT 5 (h). The repeated calculation for t = 
1000 (h) took about 1500 s with a conventional desktop PC. 
After the calculation was finished, the calculated final crystal 
size was almost the same as that we actually obtained 
(±0.02%). Therefore, we concluded that the program ran with 
satisfactory precision for the comparison of crystal growth in 
space and on the ground. 

3. Lysozyme crystal growth in THM and TDM 

3.1. Parameter calculations based on experiment results 

To use some realistic parameters, we referenced lysozyme 
crystallization using the batch method. Purified lysozyme 
(20 mg ml" 1 ) was crystallized using 0.7 M sodium chloride as a 
precipitant in 50 mM sodium acetate pH 4.5. The final size of 
the crystal [i?(oo)] and the final protein concentration in the 
solution (Ce) were measured, and the kinetic constant of the 
protein molecule (/?) and the diffusion coefficient of the 
protein molecule (D) were estimated as shown in Table 1 
(Tanaka et at, 2012). 

3.2. Lysozyme crystal growth in THM and TDM 

The time course of lysozyme crystal growth is shown in 
Fig. 2(a). The solid line is for THM obtained from the repe- 
tition of calculation (8) and the dotted line is for TDM 
obtained from the repetition of calculation of equations (16), 
(17), (20), (21), (26) and (27). This figure indicated that the 
crystal growth in TDM is a little slower than in THM. This 
might suggest that, in the case of THM, the fast uptake of the 
protein molecules into the crystal was accelerated by the 
transportation of protein molecules toward the crystal by 
convection flow. However, in the case of TDM, the PDZ 
deprived the protein molecules from around the growing 
crystal and only the diffusive flow could deliver protein 
molecules to the crystal, so it kept the protein concentration 
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Calculated results of the lysozyme crystal growth in the salt condition in the homogeneous (solid line) and diffusive (dotted line) model in the transient 
state, (a) Time course of the crystal growth, (£>) Time course of the protein concentration on the crystal surface, (c) Supersaturation level around the 
growing crystal. R(t), R: radius of the crystal (mm); C(t): concentration of the protein on the surface of the crystal (mg mP 1 ); a(R): supersaturation level; 
V. time after the crystal starts growing. 



around the surface of the crystal lower than in THM and 
caused slower crystal growth. 

3.3. Average protein supersaturation level 

The degree of supersaturation on the surface of the crystal 
while it is growing is defined as 

a(R) = [C(R) - Ce]/Ce, (29) 

where C(R) is the concentration of the protein on the surface 
of the crystal when the crystal radius is R. It is said that the 
quality of the crystal, indexed by the X-ray diffraction reso- 
lution, mosaicity and/or ^ mer ge; is better if it grows at the lower 
supersaturation (McPherson, 1999) although dependency of 
those on <r(R) has not been verified yet. To know the super- 
saturation level on the surface of the crystal while it is 
growing, the time course of the protein concentration on the 
crystal surface was calculated as shown in Fig. 2(b). The solid 
line is for THM obtained from equation (6) and the dotted line 
is for TDM obtained from equation (27). Figs. 2(a) and 2(b) 
were combined to create Fig. 2(c) using equation (29). The 
solid line and dotted line are for THM and TDM, respectively. 
It was found that the supersaturation level was high in the 
centre of the crystal and gradually became lower toward the 
surface of the crystal. Over the full range of crystal growth, the 
supersaturation level was slightly lower in TDM than in THM. 



At the end of crystal growth, the supersaturation level was the 
same as the protein solubility both in THM and TDM. 

Since X-ray diffraction intensity depends on the volume of 
the protein molecules in a crystal (McPherson, 1999), the 
quality of an X-ray diffraction image may depend on the 
integrated diffraction images obtained from a certain volume 
of the crystal which was grown in changeable supersaturation. 
Therefore, as a quantitative index of supersaturation, the 
average supersaturation level (ASS) was defined as the inte- 
grated a as the crystal grew, averaged by the total volume, 

J?(oo) 

ASS = (1/V) f a(R) dV, (30) 

R = 0 

where (j(R) is the degree of supersaturation at the surface of 
the crystal when the crystal radius is R. 

To compare the ASS within one crystal, the crystal was 
divided into three sections by volume from the centre of the 
crystal. The radius of the inner third of the crystal is 69.3% of 
R(oo), the middle 87.4%, and the outer 100%. ASS in THM 
and TDM are shown in Fig. 3(a). The ASS for a full sphere in 
THM was 1.73 and in TDM was 1.50 which was about 87% of 
that of the THM. In each third of the crystal volume, TDM 
was lower than THM. These findings are consistent with the 
estimations obtained by former analyses using the steady state 
model (Tanaka et al, 2004; Inaka et ah, 2012). 
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0.15 



s>o.io 



qO.05 



0.00 




3.4. Average impurity concentration 

As for impurity, it is not realistic to 
qualify and quantify the impurity 
molecules in a protein solution. There- 
fore, we assumed a contaminating 
impurity, whose initial concentration, 
Ci(Q), was 1% of that of the protein, 
C(0), with a diffusion coefficient, Di, the 
same (0.36 mm 2 h _1 ) as that of the 
protein, D, and a kinetic constant, f)i, 
ten times (3.4 mm h _1 ) as large as that 
of the protein fi. Di was equal to D 
because we assumed that the molecular 
weight of the impurity and the protein 
were the same. A fti ten times larger 
meant that the impurity molecule had a 
ten times higher affinity to the crystal 

than the protein molecule. Actually, when acetylated lysozyme 
was an impurity, it was reported that the impurity attached to 
the crystal with an affinity several multiples of ten times higher 
than the lysozyme protein (Thomas & Chernov, 2001). 

The time course of the concentration of the impurity around 
the surface of the crystal is shown in Fig. 4(a). The solid line is 
for THM obtained from equation (10) and the dotted line is 
for TDM from equations (16), (17), (24), (25), (26) and (28). 

The concentration of the impurity attached to the crystal 
with a radius R from the centre of the crystal [Q' cry st(^)] can be 
described as 



120 



80 120 

t(h) 

(a) 



160 



200 




0.2 0.3 
R(mm) 

(b) 



G cryst (i?) = Xi(t)/dV(t) 

= fii Ci(l, f)/{/3[C(l, t) - Ce]/n} 
= piCi(l,t)At/AR(t), 



where Xi(t) is the weight of the impurity attached to the 
crystal surface in a unit of time, and C(t) and Ci(t) are 
concentrations of the protein and the impurity on the surface 
of the crystal at time t. Generally, if many impurity molecules 
attach to a crystal, a highly disordered crystal may grow. 
Therefore, to determine the impurity concentration inside the 
crystal, Figs. 2(a) and 4(a) were combined and equation (30) 
was applied to create Fig. 4(b). The solid and dotted lines are 
for THM and TDM, respectively. As shown in Fig. 4(b), if 
the radius of the crystal was between 0 and 0.25 mm, the 
concentration of the impurity on the surface of the crystal in 
TDM was much lower than that of THM. Then, if the radius 
was larger than 0.25 mm, the impurity concentration in TDM 
became higher than that of THM, but much lower than that at 
the beginning of the crystal growth, and the concentration of 
impurity around the crystal became zero at the end of crystal 
growth. 

Similar to the ASS, the average impurity concentration 
(AIC) was defined as 



Figure 4 

Calculation results of the impurity concentration in salt conditions in the homogeneous (solid) and 
diffusive (dotted) model in the transient state, (a) Time course of the impurity concentration on the 
crystal surface, (b) Impurity concentration around the growing crystal. R: radius of the crystal (mm); 
Ci(i): concentration of the impurity on the surface of the crystal (mg mf 1 ); C;' cryst (R): concentration 
of the impurity in the crystal at R from the centre of the crystal; t: time after the crystal starts 
growing. 



As in the ASS, the AIC was calculated for a full sphere crystal 
and in each of the three sections as in Fig. 3(b). The AIC 
of TDM and THM were the same for a full sphere crystal 
because all of the impurity was finally adsorbed into the 
crystal. In Fig. 3(b), the AIC was high in the inner third of the 
crystal and fell as the crystal grew in both TDM and THM. In 
the inner third of the crystal the AIC was lower in TDM than 
in THM, but in the middle section of the crystal the AIC was 
higher in TDM than in THM. In the outer section of the 
crystal, the impurity concentration was almost zero in THM 
and TDM. These might suggest that, in TDM, the IDZ was 
formed around the crystal in the diffusive field and deprived 
the impurity molecules around the crystal surface, but, in 
THM, the fast attachment of the impurity molecules to the 
crystal lowered the concentration of the impurity more 
quickly. These findings could not be elucidated by the steady 
state model performed previously. 



(31) 



R(oo) 

aic = (i/v) / a, 

R = 0 



cryst 



(R) dv. 



(32) 



4. Conclusion 

From the results of the model calculation, it became clear 
quantitatively that all the sections of the crystal grown in THM 
and TDM are surrounded by different supersaturation levels 
of protein and different concentrations of impurity when they 
grow. Therefore, for X-ray diffraction experiments, we 
recommend considering how the quality difference within a 
whole crystal, caused by the difference in protein super- 
saturation and impurity concentration levels during crystal 
growth, may affect the quality of the X-ray diffraction 
patterns. 

As for protein supersaturation, it became clear that, in a 
diffusive model such as in microgravity, protein crystals grow 
in a lower supersaturation level of protein than in a homo- 
geneous field. This may be due to the effects of the PDZ. 

Regarding impurity, its concentration is higher in THM 
than in TDM when the inner section of the crystal grows. 
In the middle section, although the impurity concentration is 
higher in TDM than in THM, it is much lower than in the inner 
section both in THM and TDM. This may be due to the effects 
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of IDZ and the faster uptake of the impurity molecules into 
the crystal. 

It was reported that the growth of the middle and outer 
sections of a crystal are influenced by the molecular order of 
the inner section of the crystal as shown with X-ray tomo- 
graphy (Sawaura et al, 2011). Therefore, we speculate that the 
outer section of crystals grown in TDM may have better 
quality because the supersaturation ratios and impurity 
concentrations are lower in the inner section of the crystals 
grown in TDM than those of the crystals grown in THM. This 
may explain the reason why some better quality crystals grew 
in microgravity (Tanaka et al, 2007, 2011; Takahashi et al, 
2010). Further examination is required to know the influence 
of the ASS and AIC on the X-ray diffraction patterns quan- 
titatively. 

In this study, we have fixed fi as a constant. Actually the 
increase of the impurity concentration reduces the velocity of 
the crystal growth (Nakada et al, 1999). Thus, for the next 
step, we will proceed to analyze numerically how the diffusive 
model of protein crystal growth is affected by the increase of 
the initial amount of the impurity. 

Here we have introduced a practical mathematical model 
which can be calculated with an ordinary personal computer. 
Using this transient model, further numerical analyses of 
higher viscosity of the crystallization solution (lower D) and of 
higher homogeneity of the protein sample (higher f3) will be 
examined for the verification of the effectiveness of these 
values on the improvement of crystal quality in microgravity; 
and, in the near future, we will use this model on the results 
obtained in the JAXA PCG crystallization experiments to 
further understand the effects of microgravity. 

We are thankful to Professor Atsushi Nakagawa of Osaka 
University, Professor Yoshiki Higuchi of University of Hyogo, 
and Professor Sam-Yong Park of Yokohama City University in 



the International Space Station Applied Research Partnership 
Program since 2004. 
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